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ABSTRACT 

We consider the further development of the formalism of the estimates of higher-order pertur- 
bative corrections in the Euclidean region, which is based on the application of the scheme- 
invariant methods, namely the principle of minimal sensitivity and the effective charges ap- 
proach. We present the estimates of the order O^a^) QCD corrections to the Euclidean quan- 
tities: the e"'"e~-annihilation D-function and the deep inelastic scattering sum rules, namely 
the non-polarized and polarized Bjorken sum rules and to the Gross-Llewellyn Smith sum 
rule. The results for the D-function are further applied to estimate the O(a^) QCD correc- 
tions to the Minkowskian quantities R{s) = atot{^'^^~ hadrons)/(T(e"'"e" /^^A^") and 
Rr = r(r ^ u.r + hadrons)/r(r —>■ UrVec). The problem of the fixation of the uncertainties due 
to the 0(0;^) corrections to the considered quantities is also discussed. 
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I. The main aim of this work, which is based on the considerations of Refs. |1|, 0, is 
to obtain the estimates of the higher-order QCD corrections to R{s), and deep-inelastic 
scattering sum rules using the improvement formula P] of re-expansion of the expressions for 
the physical quantities obtained within the principle of minimal sensitivity (PMS) and 
the effective charges (ECH) approach , which is equivalent a posteriori to scheme-invariant 
perturbation theory 0. 

Of course, this approach to the estimation of the values of the uncalculated higher-order 
terms can not be considered as an alternative to direct analytical or numerical calculations. 
However, we hope that this method can give an idea of the possible values of the higher-order 
terms. This hope is based on the fact that the application of this method for the estimate of 
the four-loop corrections to the expression for {g — 2)e gave results which turned out to be in 
surprisingly good agreement with the latest results of the direct numerical calculations of Ref. 
[0 |. The second argument in favour of this procedure came from its successful application 
for the analysis of the Drell-Yan cross section at the 0(Q;^)-level 0. Moreover, as will be 
demonstrated in our work, the re-expansion formalism of Ref. ^ is also working quite well in 
QCD at the NNLO for at least three independent quantities, namely for the e"'"e~-annihilation 
D-function and the non-polarized and polarized Bjorken sum rules. Therefore, we prefer to 
consider these five facts not as a surprising accidental coincidence, but as an argument in favour 
of more detailed studies of the intrinsic features of this approach and its various applications. 

In this work we further develop the re-expansion formalism of Ref. Q, deriving new terms 
in the corresponding improvement formula. The previously-known terms are used to obtain 
the estimates of the next-to-next-to-next-to-leading order {N^LO) QCD corrections to the D- 
function, the non-polarized Bjorken sum rule and polarized Bjorken sum rule, which is closely 
related by the structure of the corresponding perturbative series to the Gross-Llewellyn Smith 
sum rule. We will use the results obtained for the D-function to estimate the effects of the 
N^LO corrections to the perturbative series for the Minkowskian quantities R{s) and Rr by 
adding the explicitly calculable terms, previously discussed in Refs. [p!0|, supplementing thus 



the related considerations of Refs. [|T2| by the additional input information. The derived 
new term in the improvement formula is applied to touch on the problem of fixing the values 
of the 0(a;g )-corrections to the analysed quantities. 

II. Consider first the order O(a^) approximation of a Euclidean physical quantity 

N-l 

DN = doa{l+ ^ dia') (1) 

i=l 

with a = Ofs/TT being the solution of the corresponding renormalization group equation for the 
/3-function which is defined as 

da 

/x^— = /3ia) = -/Soa'il + ^ c.a') . (2) 
(^f^ i=i 



detailed re-consideration of the analysis of Ref. and its generalization to five- loop order will be 
presented in a more detailed work m. 
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In the process of the concrete calculations of the coefficients di,i > 1 and Ci,i > 2, the MS 
scheme is commonly used. However, this scheme is not the unique prescription for fixing the 
RS ambiguities (for the recent discussions see e.g., Ref. ||T2[| ). 



The PMS 1^ and ECH prescriptions stand out from various methods of treating scheme- 
dependence ambiguities. Indeed, they are based on the conceptions of the scheme-invariant 
quantities, which are defined as the combinations of the scheme-dependent coefficients in Eqs. 
(1) and (2). Both these methods pretend to be the "optimal" prescriptions, in the sense 
that they might provide better convergence of the corresponding approximations in the non- 
asymptotic regime, and thus allow an estimation of the uncertainties of the perturbative series in 
the definite order of perturbation theory. Therefore, applying these "optimal" methods, one can 
try to estimate the effects of the order 0(a^"'"^)-corrections starting from the approximations 
D°^^{aopt) calculated in a certain "optimal" approach 0], 0, [1^ . 



Let us following the considerations of Ref. and re-expand D°^^{aopt) in terms of the 
coupling constant a of the particular scheme 

DtM = D^ia) + a^+i (3) 

where 

5Dt = nNid,,c,) + n^id:i^\cf) (4) 

are the numbers which simulate the coefficients of the order 0(a^^^ )-corrections to the physical 
quantity, calculated in the particular initial scheme, say the MS'-scheme. 

The explicit form of the coefficients fli can be obtained following the considerations of Ref. 
0. We present here the final already known expressions 0: 

^2 = dodi{ci + di), (5) 

= dodi{c2 - ^cirfi - 2dl + 8(^2) . (6) 
and the new term evaluated by us: 

= — (Scgdi + 02^2 - 4:C2dj + 201^1^2 - Cirfg + Ud'^ - 28djd2 + 5dl + Udids) (7) 
3 

which reproduces the renormalization-group controllable logarithmic terms at the five-loop level 

It should be stressed that in the ECH approach df^^ = for alH > 2. Therefore one gets 
the following expressions for the NNLO and N^LO corrections in Eq. (3): 

5Df'=^^ = fi2(^^l,Cl) (8) 

5Di^" = n,{d,,d2,Cr,C2) (9) 
(5Df = fi4(cii, 4, C^3, Cl, C2, Cs) (10) 
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In order to find similar corrections to Eq. (3) in the N-th order of perturbation theory starting 



from the PMS approach Q], it is necessary to use the relations obtained in Ref. |]T5| between 
the coefficients rf^^ and cf*^"^ {i > 1) in the expression for the order 0{a^j^^.^g) approximation 
D^'^^iapMs) of the physical quantity under consideration. The corresponding corrections have 
the following form: 

5Df = 5Df + ^ (11) 

5D^^' = SDi'''' . (12) 

Notice the identical coincidence of the N^LO corrections obtained starting from both the PMS 
and ECH approaches. A similar observation was made in Ref. using different (but related) 
considerations. 

In the fourth order of perturbation theory the additional contribution to 6D^^''^ has more 
complicated structure. Indeed, the expression for Q^i^df^^ ,cf''^^^) in Eq. (4) reads: 

n,{dr', cr') = y [^c,cr^ - l(cr^)^ - |c?cr*^^ + ^cf] (13) 



where 



and 



cf ^5 = ^(ci2 + C2 - ci? - c,d, + ^cl) + 0{a) (14) 
8 3o 



cf^« = 4(^3 + - C2d, - 3did2 + 2d\) + ]-Ci{d2 + C2 + Ml - cidi + -^cj) + 0(a). (15) 

/ / iUo 

The expressions for Eqs. (14)- (16) are pure numbers, which do not depend on the choice of 
the initial scheme. Note, that we have checked that in the case of the consideration of the 
perturbative series for {g — 2)e and {g — 2)^ the numerical values of fl^^df^^, cf^^) are small 
and thus the a posteriori approximate equivalence of the ECH and PMS approaches is preserved 
for the quantities under consideration at this level also We think that this feature is also 
true in QCD. 

III. Consider now the familiar characteristic of the e~^e~ — > 7 — > hadrons process, namely 
the D-function defined in the Euclidean region: 



00 



R{s 



Its perturbative expansion has the following form: 

D(g') = 3^Q}[l + a + J2dia'+'] + {J:Qff[d2a' + 0{a^)] (17) 

i>l 

where Qf are the quark charges, and the structure proportional to (SQ/)^ comes from the 
light-by-light-type diagrams. The coefficients di and d2, d2 were calculated in the MS'-scheme 
in Refs. [|l6l and |18| respectively. Following the proposals of Ref. [|l9l, we will treat the 



light-by-light-type term in Eq. (18) separately from the "main" structure of the D-function, 
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which is proportional to the quark-parton expression D^^IQ"^) = SSQj. In fact, one can hardly 
expect that it is possible to predict higher-order coefficients di,i>3 of the second structure in 
Eq. (17) using the only explicitly- known term d2- Therefore we will neglect the light-by-light- 
type structure as a whole in all our further considerations. This approximation is supported 
by the relatively tiny contribution of the second structure of Eq. (17) to the final NNLO 
correction to the D-function. 

The next important ingredient of our analysis is the QCD /5-function (2), which is known 
in the MS-like schemes at the NNLO level |2|]. 



Using now the perturbative expression for the D-function, one can obtain the perturbative 
expression for R{s), namely 

R{s) = 3SgJ[l + a, + 5^ r,ai+'] + {J^Qffihal + ...] (18) 

i>l 

where = Os/vr, ri = di, r2 = d2 — h'^Pq/S, f2 = and = d^ — n'^PQ^di + |ci) . The 
corresponding vr^-terms come from the analytic continuation of the Euclidean result for the 
D-function to the physical region. The effects of the higher-order vr^-terms were discussed in 
detail in Ref. [|1^]. For example, the corresponding expression for the r4-term reads : 

7 1 TT^ 

r4 = ^4 - vr2/?2(2rf2 + -cirfi + -cl + C2) + y/5o" (19) 



The perturbative expression for Rr is defined as 

f-M? ds 

(1 - s/M;y (1 + 2s/M;)R[s) ~ 3[1 + + 

i>l 



Rr = 2j^ ^ (1 - s/MlY (1 + 2s/Ml)R{s) ^ 3[1 + a. + E rja^^^'] (20) 



where a^- = as{M^)/TT and R{s) is R{s) with with / = 3, (SQ/)^ = 0,3T,Q'j substituted for 
3E I Vff> p and | P + | Ks 1. 

It was shown in Ref. |Tl| that it is convenient to express the coefficients of the series (24) 



through those ones of the series (18) for the D-function in the following form : 

r[ = df^{f = 3)+g, 

rl = rfr(/ = 3)+^72 

rl = df'{f = 3)+gs (21) 

where in our notations 

gi = -Poh = 3.563 

92 = -[2d, + c,]/3oh + f3ll2 = 19.99 

93 = -[34 + 2ciiCi + C2]/3o/i + [3cii + ^01^/2-/33/3 = 78.00 (22) 
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and Ik are defined and calculated in Ref. [^. Their analytical expressions read: /i=-19/12, 
/2=265/72-7rV3 and /3=-3355/288+ 197rVl2. One of the pleasant features of Eqs.(^ is that 



they are absorbing all effects of the analytical continuation. 

Following the lines of Ref. |TI| we derive the corresponding expression for the coefficient rl: 

rl = df^if = 3)+g, (23) 

where 

3 

§4 = -[Ads + 3d2Ci + 2diC2 + C3]l3oh + [Qd2 + 7cidi + -cI + 3c2]I3qI2 
13 

-[Ad, + Y^iKh + P^h = 3.562c3(/ = 3) + 14.2474(/ = 3) - 466.818 (24) 

and /4 =41041/864-2657rV36+7r75 ^ -5.668. 

In order to estimate the values of the order O(a^), 0{a'^) and O(a^) corrections to R{s) 
and Rr, we will apply Eqs. (5) - (10) in the Euclidean region to the perturbative series for the 
Z)-function and then obtain the estimates we are interested in using Eqs. (18), (19), (21)-(24). 
This is the new ingredient of this analysis, which distinguishes it from the related ones of Refs. 

i, El- 



We now recall the perturbative expression for the non-polarized Bjorken deep-inelastic scat- 
tering sum rule 

BjnSR = C F^'^^ix, Q^)dx = 1 - -a(l + V d,a') (25) 
Jo 3 XI 



where the coefficients di and d2 are known in the MS scheme from the results of calculations 
of Ref. pri and Ref. P5[ . The expression for the polarized Bjorken sum rule BjpSR has the 



following form: 

^ ep—eni /W2\ j ^ I 9^ 



BjpSR = [ (7r"""(x, Q^)dx = - I — I [1 - a(l + V d,a')] (26) 
^0 3 gv 



where the coefficients di and d2 were explicitly calculated in the MS scheme in Refs. and 
[El] respectively. 



It is also worth emphasizing that, in spite of the identical coincidence of the NLO correction 
to the Gross-Llewellyn Smith sum rule GLSSR = (1/2) F^'^^^^{x,Q'^)dx with the ones for 
the BjpSR the corresponding NNLO correction differs from the one for the BjpSR by 



the contributions of the light-by-light-type terms typical of the GLSSR [^]. Since these light- 
by-light-type terms appear for the ffist time at the NNLO, it is impossible to predict the 
value of the light-by-light-type contribution at the N^LO level using the corresponding NNLO 
terms as the input information. However, noticing that at the NNLO level the corresponding 



light-by-light-type contributions are small |24], we will assume that the similar contributions 



are small at the N^LO level also. Only after this assumption can our estimates of the NNLO 
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and N^LO corrections to the BjpSR be considered also as the estimates of the corresponding 
corrections in the perturbative series for the GLSSR . 

IV. The estimates of the coefficients of the order O(a^), 0{a'^) and O(a^) QCD correc- 
tions to the £)-function, R{s), BjnSR and BjpSR/GLSSR obtained from the ECH-improved 
expressions are presented in Tables 1-4 respectively. Due to the complicated /-dependence of 
the coefficients Q2,^3 , ^4 in Eqs. (5)-(7), we are unable to predict the explicit /-dependence 
of the corresponding coefficients in the form respected by perturbation theory. The results 
are presented for the fixed number of quark flavours 1 < / < 6 and are normalized to the 
M5'-scheme. 



/ 






df^ 


df^ - czdi 


1 


14.11 


7.54 


75 


474 


2 


10.16 


6.57 


50 


261 


3 


6.37 


5.61 


27.5 


111 


4 


2.76 


4.68 


8 


23 


5 


-0.69 


3.77 


-8 


-15 


6 


-3.96 


2.88 


-21 


-1.8 



Table 1 : The results of estimates of the NNLO, N^LO and N'^LO corrections in the series 
for the D-functions. 



/ 






^est 


* - Cgdi 


1 


-7.84 


-14.41 


-166 


-1748 


2 


-9.04 


-12.63 


-147 


-1156 


3 


-10.27 


-11.03 


-128 


-669 


4 


-11.52 


-9.58 


-112 


-263 


5 


-12.76 


-8.29 


-97 


64 


6 


-14.01 


-7.17 


-83 


334 



Table 2: The results of estimates of the NNLO, N^LO and N'^LO corrections in the series for 
R{s). 



6 



f 

J 






df* 


df * — C'^d^ 


1 


44.97 


39.62 


424 


4177 


2 


36.19 


33.28 


303 


2619 


3 


27.89 


27.37 


200 


1481 


4 


20.07 


21.91 


114 


671 


5 


12.72 


16.91 


44 


137 


6 


5.85 


12.39 


-10 


-148 



Table 3: The results of estimates of the NNLO, N^LO and N'^LO corrections in the series for 
BjnSR. 



/ 


df 


df* 


df* 


df * - Cadi 


1 


34.01 


27.25 


290 


2557 


2 


26.93 


23.11 


203 


1572 


3 


20.21 


19.22 


130 


854 


4 


13.84 


15.57 


68 


342 


5 


7.83 


12.19 


18 


27 


6 


2.17 


9.08 


-22 


-135 



Table 4: The results of estimates of the NNLO, N^LO and N^LO corrections in the series for 
BjpSR and GLSSR. 

The presented NNLO estimates are in reasonable agreement with the results of the explicit 
calculations. The best agreement is achieved for / = 3 numbers of flavours. It can be shown 
that the difference with the results obtained staring from the PMS improved expression is rather 
small Big. 

In order to address the problem of fixing the values of the O(a^) QCD corrections to the 
considered Euclidean quantities, Eqs. (7), (10) with explicitly calculated NLO and NNLO 
coefficients di,d2,Ci,C2 were applied. The values of the coefficients ds were fixed by using the 
estimates of the N^LO coefficients the determined from Eqs. (6), (9). To estimate the next- 
to-next-to-next-to-next-to-leading order (N'^LO) coefficient of R{s), which is related to the 
N^LO coefficient of the D-function, the explicitly calculated terms in Eq.(19) were taken 
into account. However, the expressions for Q4 in Eq. (7) depend also on the four-loop coefficient 
C3 of the QCD /3-function which is unknown at present. Therefore, in Tables 1-4 we presented 
the estimates of the combinations d4 — diCs and — ric^. Another existing uncertainty in the 
value of N^LO coefficients da is fixed by the assumption that the real values of these coefficients 
do not significantly differ from the N^LO estimates obtained by us. However, in the case of 
the D-function with / = 3 numbers of flavours we present the more detailed expression of d4: 

df *(/ = 3) = 1.64c3(/ = 3) + 5.97d3(/ = 3) - 52.66. (27) 

In order to obtain the corresponding estimate from Table 1, Eq.(27) should be supplemented by 
the A^^LO-estimate df *(/ = 3) = 27.5 flrst obtained in Ref.||. This result is in good agreement 
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with the "geometric progression" assumption d^^f = 3) ~ dl{f = 3)/di{f = 3) ~ 25 of Ref. ||TT 



Our result was recently supported ||2^ by the phenomenological analysis of the ALEPH data 
for Rt-. 

Amongst other most interesting results of our analysis are the estimates of the N^LO and 
N^LO corrections to R{s) for / = 5 numbers of flavours. Taking as{Mz) ~ 0.12, we get the 
estimate of the corresponding N^LO contribution to r{Z^ hadrons) 



i5Tzo)N^LO ~ -97{a{Mz)Y 



-2 X 10" 



(28) 



The comparison of the results of Table 1 with those of Table 2 demonstrate that the vr^-effects 
give dominating contributions to the coefficients of R{s) and F^o. Indeed, at the A^^LO-level 
they are also larger than the estimated coefficient d^ of the Euclidean D-function. Assuming 
that cs{f = 5) = C2(/ = 5)^/ci(/ = 5) ~ 1.715 we arrive at the following less substantiated 
estimate of the corresponding N'^LO contributions, namely 



N-^LO 



64:{a{Mz) f ^ 5 X 10 



-6 



The way of fixation of the value of C3 used above is not applicable for the case of / 
we expect that in this case its real value is negative. 



(29) 



6, since 



Using now the relations derived in Refs. pq , |rT| between the coefficients of the D-function 
and Rr and taking as(M^) ^ 0.36 |^ we get the following estimate of the N^LO contribution 
to Rt-: 

((5i?^)^3io ~ 105.5 ^ 1.8 X 10^2 ^ (30) 



This estimate is more precise than those presented in Refs. and 0], namely 6Rr = ±130 



and 6Rr = (78±25)a^ |TT]], and is smaller than the result of applying the Fade resummation 
technique directly to Rr, namely 6Rr = 133a^ ||27[| . 



Generalizing the considerations of Ref. |TT| to the order 0(a^)-level and taking c^lf = 3) ~ 
C2(/ = 3)^/ci(/ = 3) ~ 11 we get the estimate of the N^LO corrections to Rr 0). This estimate 
reads 01: 



iSRr] 



N^LO 



Mai 



1.8 X 10' 



(31) 



V. Let us emphasize that the results for R{s) and Rr were obtained by us from those of the 
D-function after taking into account explicitly calculable terms of the analytic continuation to 
the Minkowskian region, discussed in detail in Ref. In principle, one can try to study the 
application of the procedure used by us directly to R{s) and Rr in the Minkowskian region. We 
checked that in the case of applications of Eqs. (5)- (7) with the substitution of di = ri or di = rj 
one can reproduce the same values of the Eucledian coefficients di but the structure of the high 
order vr^-terms will not agree with the explicitly known results. A similar problem will also 
definitely arise in the case of more rigorous studies of the applicability of the Fade resummation 
methods for the Minkowskian quantities R{s) and Rr (for the studies existing at present see 
Refs. [0). We think that these analyses should be supplemented by the construction of the 
Fade approximants directly to the D-function. In our case it is possible to reproduce our results 
R{s) and Rr by means of application of a similar technique directly in the Minkowskian region 
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after modifications of Eqs. (5)-(7) in the Minkowskian region by adding concrete calculable 
TT^-dependent and scheme-independent factors. 
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